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Abstract 

A formalism for anisotropic fluid dynamics is proposed. It is designed to describe 
boost-invariant systems with anisotropic pressure. Such systems are expected to be 
produced at the early stages of relativistic heavy-ion collisions, when the timescales 
are too short to achieve equal thermalization of transverse and longitudinal degrees 
of freedom. The approach is based on the energy-momentum and entropy conser- 
vation laws, and may be regarded as a minimal extension of the boost-invariant 
standard relativistic hydrodynamics of the perfect fluid. We show how the formal- 
ism may be used to describe the isotropization of the system (the transition from the 
initial state with no longitudinal pressure to the final state with equal longitudinal 
and transverse pressure). 
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1. At present, the most successful description of early parton dynamics is 
achieved with the help of the relativistic hydrodynamics [Tf2"f5] . With the 
equation of state incorporating the phase transition and with the appropriate 
modeling of the freeze-out process, the hydrodynamic approach leads to very 
successful description of the hadron transverse-momentum spectra and the 
elliptic flow coefficient v 2 [4151Id] . We also note that with a suitable modification 
of the initial conditions, the hydrodynamic approach describes consistently the 
HBT radii [7]. 
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In spite of those clear successes, the use of the hydrodynamics is faced with 
the problem of so-called early thermalization - in order to have a successful 
description of the data, the hydrodynamic evolution (implicitly assuming the 
three-dimensional local equilibrium) should start at a very early time, well 
below 1 fm/ c. Such short values can be hardly explained within the microscopic 
calculation. 

Recently, a possible solution to the problem of early thermalization has been 
proposed [8] . With the assumptions that only transverse degrees of freedom are 
thermalized and the longitudinal dynamics is essentially the free-streaming, 
one can obtain the parton transverse-momentum spectra and V2 which agree 
well with the data [8119] . In this approach, called below the transverse hydro- 
dynamics (for general formulation see [ID]), the longitudinal pressure vanishes 
while the transverse pressure is large and leads to the formation of a substan- 
tial transverse flow, which is the main effect responsible for the good agreement 
with the data. 

One naturally expects that after some time the purely transverse hydrody- 
namic evolution is transformed into the standard hydrodynamic evolution with 
isotropic pressure. The typical way to describe such transformation would be 
to use the kinetic theory [TTJ or dissipative hydrodynamics [T2p3] . Such ap- 
proaches lead, however, to the entropy and particle production which seems to 
contradict the observed scalings of the hadron production with the number of 
the initial constituents (unless the entropy production is a small effect). An- 
other mechanism to describe a similar transition, from the initial quasithermal 
two-dimensional parton distribution to the final three-dimensionally isotropic 
parton distribution, was studied in Ref. [H]. This approach conserves the en- 
tropy and is based on the coupled Vlasov and Yang-Mills equations for the 
quark-gluon plasma. The full isotropization of the system is in this case an 
effect of bending of parton trajectories in strong color fields. 

The aim of this paper is to propose the extension of the standard boost- 
invariant hydrodynamics, which would be suitable for the description of sys- 
tems with anisotropic pressure. In the special cases, our approach is reduced 
to the transverse hydrodynamics or standard hydrodynamics. It may be also 
used to describe effectively the process of full isotropization of the pressure. 

The proposed formalism is based on the energy, momentum, and entropy con- 
servation laws. As we have mentioned above, the use of the entropy conser- 
vation is suggested by various modeling of the RHIC data. For example, the 
PHOBOS data [15] shows that the numbers of produced hadrons per wounded 
nucleon in central d+Au and Au+Au collisions are very much similar (dif- 
ferences smaller than 30%). This indicates that equilibration/isotropization 
effects do not produce a large amount of the entropy. In addition, the re- 
cently observed scalings of the hadron production with the number of so-called 
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wounded constituents [16] leave little room for extra particle production dur- 
ing the evolution of the system. 

2. Our starting point is the following form of the energy- momentum tensor 
T» v = (e + P T ) U>*U U - P T gT - (P T ~ P L )V^V V ', (1) 



where e is the energy density, Pt and Pl are the transverse and longitudinal 
pressure, and U^ 1 is the four-velocity of the fluid satisfying the normalization 
condition U^U^ = 1. In the isotropic case, the pressures Pt and Pl are equal, 
Pt = Pl = P, and the energy-momentum tensor takes the standard form. For 
the anisotropic fluid, the last term in ([T]) is different from zero. The fourvector 
defines the direction of the longitudinal pressure. It is spacelike, orthogonal 
to £/ M , U^Vu = 0, and normalized by the condition V^V^ = —1. In the local 



rest frame of the fluid element we have = (1, 0, 0, 0) and ■ 
hence the energy-momentum tensor ([I]) takes the expected form 



(0,0,0,1), 



e ^ 
P T 
P r 
\0 P L ) 



(2) 



For the boost-invariant systems the fluid four-velocity has the structure 
£/ M = (vP cosh 7] } u X: Uyj vP sinh T^j , (3) 



where = (u°, u x , u y , 0) is the fluid four-velocity at z = (for the vanishing 
longitudinal coordinate), and r] = 1/2 \n(t+z)/ (t—z) is the spacetime rapidity. 
With the normalization condition Uq — u x — u 2 = 1 we automatically have 
U 2 = 1. The four- vector satisfying the appropriate normalization and 
orthogonality conditions has the form 

V 1 = (sinh T), 0,0, cosh rj) . (4) 



In this paper we restrict our considerations to the case of massless particles, 
where = and e = 2Pt + Pl- As in the standard hydrodynamics, the 



Throughout the paper we use the natural units with c = h = 1. The metric tensor 
= diag(l, -1,-1,-1). 
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evolution equations are obtained from the energy-momentum conservation 
law 



d,T^ = 0. (5) 

The projection of Eq. (jSJ) on the four-velocity U v yields 

De + (e + P T )d^ - AU v V t d,y v = 0. (6) 

where we have introduced the short-hand notation for the total time derivative, 
D = U^d^ and the difference of the pressures, A = P T — P L . In addition to the 
energy-momentum conservation law (J5J) we demand that there is a conserved 
entropy current characterizing the system. We write it in the form d^aU^) = 
or equivalently as 

Da + ad^U^ = 0. (7) 

where a is the entropy density. Moreover, with the definition (j3J) one finds 
V^V = ff' hir, (8) 



where r = y/t 2 — z 2 is the longitudinal proper time. Eqs. (JTj) and (jSJ) allow us 
to rewrite Eq. (jSJ) in the form 

De= (£±M Da + ^ DT . (9) 

a r 



Eq. (jUJ) indicates that in the general case the energy density may be considered 
as a function of the two variables, e = e(a, r), hence e(r, x, y) = e [cr(r, x, y),r]. 
We emphasize that this is a novel feature of our approach, which distinguishes 
it from the standard hydrodynamics where the energy density depends only 
on the entropy density, e = e(a) and e(r, x,y) — e [cr(r, x, y)]. 

3. The functional dependence e = e(a, r) plays a role of the generalized equa- 
tion of state for our system. Eq. (jUj) is satisfied if the two conditions hold, 



Multiplying Eqs. ffTUl) by a and r, respectively, and dividing both of them by 
e we obtain a simple set of equations 

(de'\ 4 A' (de'\ 
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where e' = ln(e/eo)) A' = A/e, cr' = ln(cr/cr ), and r' = ln(r/r ) with e ,o"o,' r o 
being arbitrary constants UJ. The thermodynamic consistency requires that 
de' is the total derivative, hence the mixed second derivatives of the function 
e'{a' , t') are equal. This condition implies that A' must be a function of the 
single variable, namely A' = A'(o"' + 3r'). By the direct integration of Eq. (jlip 
we find that the energy density function must be of the form 



a \ 4 / 3 
e \ — ) R(x) 



x 



IT 

0"o V^o 



(12) 



The function R(x) is related to the function A'(x) by the equation 

In x 



R(x) = exp 



1 



A'(y)dy 



(13) 



In the similar way we may express the transverse and longitudinal pressure, 
namely 



. ° \ 4/3 
Pt = s ( — 
-cr / 



R 



+ xR' 



a \ 4 /3 
Pl = s (-) 



R 



-2xR' 



(14) 



where R' = dR/dx. We note that Eqs. (|T2|) and (fll"]) satisfy explicitly the 
condition £ = 2P^ + P^. 

Equation fflBl indicates that R depends on the space-time evolution of the 
pressure anisotropy. As a consequence, its particular form may be obtained 
only from the microscopic calculation. At the hydrodynamic level of descrip- 
tion, which is the subject of our analysis, the form of the function R should 
be assumed. Nevertheless, from the constraints Pt > Pl and Pl > 0, we ob- 
tain the two conditions: R'(x) > and R(x) > 6xR'(x), hence R(x) cannot 
grow faster than x 1//6 . Moreover, we expect that R saturates for large values of 
time, when the transition to the standard hydrodynamics takes place. Below 
we shall refer to R as to the pressure relaxation function. With a given form of 
R, the generalized equation of state (fl2"|) and also the dynamics of the system 
are well defined. 



We define the temperature in the standard way as the derivative of the energy 
density with respect to the entropy density, 




(15) 



2 Their values may be arranged to impose the appropriate initial conditions. 



5 



which together with Eq. (Q leads to the following generalized thermodynamic 
identities 



Tcr = e + P T , De = TDa + AD In r, DP T = oDT — AD In r. (16) 



One can notice that for the isotropic system, where A = and Pt = Pl = P , 
the above equations are reduced to the standard thermodynamic identities 
valid for the systems with zero baryon density. 

4. Starting from the energy-momentum conservation laws (JSJ) and using Eqs. 
(1161) we obtain the final form of the hydrodynamic equations for the anisotropic 
perfect fluid: 



(TU V ) = d u T, 

M (<7tf")=O. (17) 

In the derivation of Eqs. ffT7|) . which are the main formal result of our paper, we 
have used the following properties: V^d^r = 0, d^V^ = 0, and V^d^A = 0. 
They follow directly from the definition (j4j) and the assumption of boost- 
invariance. Since the projection of the first equation in (|T7|) on U v yields iden- 
tically zero and, in addition, the system is boost-invariant, Eqs. (fT7|) consist 
of 3 equations for 4 unknown functions: u x , u y , T, and a. Similarly as in the 
standard hydrodynamics, the "dynamical" Eqs. (fTT|) must be supplied with 
the "material" equation — the equation of state in the standard hydrodynam- 
ics or the generalized equation of state of the form (fT2l) in our approach. It is 
very interesting to observe that Eqs. ffTTl) have exactly the same form as the 
standard hydrodynamic equations for boost-invariant baryon-free matter. The 
only difference remains in the use of the generalized equation of state including 
the pressure relaxation function R. The latter defines the temperature T as a 
function of a and r. Equivalently, it may be used to define a as a function of 
T and r. 



5. Let us now consider a few special cases. To illustrate the main characteristic 
features of our approach with simple and analytic examples, we neglect the 
transverse expansion of the fluid. In this case, the entropy equation has always 
the solution of the form 

a[T(r),r]=a ^. (18) 



i) The isotropic case with P T = P L corresponds to the trivial choice A' = 
and gives R(x) = 1. The generalized equation of state in this case is reduced to 
the well known black-body equation e = £o( cr / cr o) 4//3 > which yields a ~ T 3 . In 
this way we recover the famous result of the Bjorken model where T ~ t -1 / 3 . 
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ii) Another interesting situation corresponds to the case with vanishing lon- 
gitudinal pressure. For Pl = one finds A' = 1/2, which in turn leads to 
the relation R(x) = x 1//6 . In this case the generalized equation of state gives 
a ~ T 2 /t [8~1TTU] . Comparing this result to Eq. (18) we conclude that the tem- 
perature becomes independent of time. This is an expected result of transverse 
hydrodynamics - without transverse expansion the transverse clusters cannot 
cool down and their temperature remains constant. 

iii) As the third case w consider the pressure relaxation function of the form 
R(x) = [xj — (xj — 1) exp[— (x — l)/(xi — l)]] 1 . Having in mind that the 
variable x grows with time, x = (<t/<To)(t/to) 3 = (r/r ) 2 , this form of the 
function R may be useful to describe a smooth isotropization transition from 
two-dimensional purely transverse hydrodynamic expansion to the isotropic 
expansion governed by the standard hydrodynamics. The parameter Xj de- 
termines here the timescale 77 for the isotropization transition, rf = XjTq. 
For a; - 1 < - 1 we have R ~ x 1//6 , hence we recover the conditions of 
transverse hydrodynamics. On the other hand, for x — 1 ^> Xj — 1 the function 
R becomes a constant, R = x 1 / 6 , and we recover the conditions of standard 
hydrodynamics, however, with a rescaled initial energy density GO- 

iv) Finally, we consider the case R(x) = x x ^d{xi — x) + x\ Q{x — xi), where 
9{x) is the step function. Similarly to the previous case, the parameter xi 
determines the timescale for the isotropization transition. However, in this 
case we deal with a sudden isotropization transition. Certainly, this kind of 
the transition may be treated only as an approximation to real processes. 

In Fig. 1 we show four examples of the function R, introduced above in the 
cases i) - iv). The solid line describes the standard hydrodynamics (this type of 
the evolution may be considered as the case with tj = r , i.e., the isotropization 
takes place right at the beginning of the evolution), the dotted line corresponds 
to the purely transverse hydrodynamics (17 = 00, the isotropization never 
happens), while the dashed and dashed-dotted lines correspond to the forms 
discussed in the cases iii) and iv) — the smooth and sudden isotropization 
with Tj = 2t . 

In Fig. 2 we show the time evolution of the energy density determined by the 
relaxation functions shown in Fig. 1 (with the same notation). The solid line 
describes the case of the standard 1+1 boost-invariant hydrodynamics where 
e ~ t~ 4 / 3 , whereas the dotted line describes the case of the transverse hydro- 

3 If the system approaches the limit of standard hydrodynamics with local thermal 
equilibrium, the combination of the constants xy 6 eoa 4//3 must be expressed by the 
number of the internal degrees of freedom. This condition leads to the constraint 
for £0 and 00, which are the initial parameters for the transverse hydrodynamics. 
The constraint may be easily fulfilled since for the transverse hydrodynamics the 
combination eo"~ 4//3 is not a constant |10| . 
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Fig. 1. Four examples of the pressure relaxation function R(x) discussed in the text 
in the cases i) - iv). 




Fig. 2. Time evolution of the energy-density in the cases i) - iv). 



dynamics where e ~ r _1 . The dashed and the dashed-dotted lines describe 
the situations with a transition from transverse to standard hydrodynamics. 
We observe that the more the evolution resembles that of transverse hydrody- 
namics, the slower is the decrease of the energy density. This is clearly caused 
by the fact that no longitudinal work is done in the case of the transverse 
hydrodynamics, hence less energy is "lost" in this case. 

In Fig. 3 we describe the temperature evolution. Again the solid line de- 
scribes the standard 1 + 1 boost-invariant hydrodynamics where T ~ r -1 / 3 , 
while the dotted curve represents the case of transverse hydrodynamics where 
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Fig. 3. Time evolution of the temperature for the cases i) — iv). 

T = const 0- Interestingly, in the case with the sudden isotropization the tem- 
perature exhibits a discontinuity. This behaviour is caused by the fact that the 
longitudinal degrees of freedom are suddenly switched on . Consequently, the 
energy and entropy must be shared among larger number of degrees of free- 
dom and the temperature decreases (note that this behavior is a consequence 
of the assumptions on the energy and entropy conservations). 



6. We close the paper with the following conclusions: 

i) Our main formal result is the set of the hydrodynamic equations (jTT]) sup- 
plied with the generalized equation of state (fT2|) . These equations determine 
the time evolution of boost-invariant perfect fluids with anisotropic pressure. 
The generalized equation of state includes the pressure relaxation function R, 
which may be used to define the way in which the system reaches complete 
isotropization. 

We note that although the examples given in this paper illustrated the 1+1 
boost-invariant expansion, Eqs. ( 1TT1) are suitable for the numerical analysis of 
the 3+1 boost-invariant systems. 

ii) The use of the proposed formalism in the phenomenological approaches to 
relativistic heavy-ion collisions is straightforward and may bring interesting re- 
sults. For example, the presence of the initial purely transverse hydrodynamic 
stage leads to the formation of a very strong flow |8j. As it has been shown 
recently in Ref. [4], the presence of the strong flow is one of the necessary 
conditions required to describe uniformly the RHIC soft hadronic data. 



4 The difference at the beginning of the two curves is caused by our normalization 
which assumes the same energy density at r = tq 
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iii) The internal consistency of the approach requires that the pressure relax- 
ation function R(x) does not grow faster than x 1//6 , and saturates for large 
arguments (in the cases when the evolution is sufficiently long and the system 
becomes isotropic). The fitting of the heavy-ion data has the potential of de- 
termining the relaxation function R, whose behavior is determined essentially 
by one parameter - the timescale for the isotropization transition. 

iv) The entropy conservation assumed in our approach suggests that our de- 
scription of the full isotropization may be an effective way of dealing with 
the effects caused by the mean-field quark-gluon plasma dynamics, based on 
the entropy-conserving collisionless kinetic equations [T7|14II18] . This issue de- 
serves further studies. Similarly, the generalizations of our approach to non 
boost-invariant cases represents an interesting and open problem. 
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